! ***************************************
finish
/clear
/filn, E333 
/title, Buckling analysis of ellipsoidal head with vertical support load
/units,si
! ********* 趨*********
Di=1000                    !Ͳ弰Բͷھ
t=10                       !Ͳ弰Բͷں
t1=18                      !ͷ
Di=1000+t                  !о
Hc=Di                      !Ͳ峤
dr=630                     !֧װֱ
ds1=150                    !ֱ
d2=108                     !ֱֹ֧
nh=30                      !ʷ
nn=20                      !ͷʷ
fyy=4.27e6                 !ֵṹ趨ôֱغ
EXX1=2e5                   !ϵĵģ
NUXY1=0.3                  !ϵĲɱ
sy=235                     !
!****************ǰ***************************
/PREP7 
et,1,181                   !趨Ԫ
R,1,t,t,t,t                !ǵԪʵ
R,2,t1,t1,t1,t1            !ǵԪʵ
mp,ex,1,EXX1               !ϵĵģ
mp,nuxy,1,NUXY1	           !ϵĲɱ

TB,BISO,1,1,2,  	   !˫Ըͬǿ
TBTEMP,0                   !¶ֵΪ0
TBDATA,,sy,0,,,,           !
!****************ģ***************************
CSWPLA,11,1,0.5,1,         !Բϵ
K,3,-Di/2,,,               !Բ˵
K,4,Di/2,,, 
L, 3,4                     !Բ
csys,0                     !ֱϵ
k,5,Di/2,-Hc               !Ͳ˵
l,4,5                      !Ͳĸ
k,1,0,0,0                  !תĹؼ
k,2,0,hc,0                 !תĹؼ
l,1,2                      !
lptn,1,3                   !߶λ
LDELE,5, 6, ,1             !ɾ߶
LDELE,4                    !ɾ߶εĹؼ
Alls                       !ȫѡ
AROTAT,2,7 , , , , ,1, 6,60, ,	!תԲͷͲĸɷͷ
wpoff,dr/2                 !ƶƽ
wprot,0,90                 !תƽ
cyl4, , ,ds1/2, , , ,-Di   !ɸԲ

aslv                       !ѡ
asel,r,loc,z,-ds1,0        !ѡ
asel,u,loc,y,0             !ȥ¶
asel,u,loc,y,Di            !ȥ϶
*get,acon,area,,num,min    !ȡи
asel,s,loc,y,0,Di          !ѡͷ
*get,acon1,area,,num,min   !ȡͷ

aptn,acon1,acon            !滥֣26
VDELE,1, , ,1              !ɾԲ

lsel,s,loc,y,Di            !ѡи϶߶
asll                       !ѡ
ADELE,all, , ,1            !ɾѡ(9)
lsel,s,loc,y,0             !ѡy=0λ߶
lsel,r,radius,,ds1/2       !ѡи¶߶
asll                       !ѡ
ADELE,all, , ,1            !ɾѡ棨10

CYL4, , ,d2/2, , , ,-Di    !ֹ֧ܸԲ
lsel,s,radius,,d2/2        !ѡ뾶Ϊd2/2߶
asll                       !ѡ
asel,u,loc,z,0,d2/2        !ȥλô0d2/2 
*get,acon,area,,num,min    !ȡиţΪ5
alls
APTN,7,acon                !ָ
VDELE,1,,,1                !ɾԲ
ADELE,10,11,,1             !ɾи

!**************************************************
lsel,s,loc,y,-hc,0         !ѡλy=-hc0֮߶
lsel,u,loc,y,-hc           !ȥλy=-hc߶
lsel,u,loc,y,0             !ȥλy=0߶ΣѡlineΪ12
LESIZE,all, , ,nh,, , , ,1 !趨Ͳ徭ʷ
lsel,s,loc,y,-hc           !ѡλy=-hc߶
lsel,a,loc,y,0             !ѡλy=0߶ΣѡlineΪ45
LESIZE,all, , ,nn,, , , ,1 !軷ʷ
alls
LESIZE,3, , ,2*nn,, , , ,1 !߶ʷ
LESIZE,19, , ,nn,, , , ,1 
LESIZE,20, , ,nn,, , , ,1 
LESIZE,16, , ,nn,, , , ,1 
LESIZE,22, , ,3,, , , ,1
LESIZE,23, , ,3,, , , ,1  
LESIZE,18, , ,2*nn,, , , ,1
LCCAT,3,5                  !ε߶Ϊһ߶

wprot,0,0,90               !תƽ棬Աʷ
ASBW,6                     !ʷֹ֧ڵ
AMESH,1                    !Ͳӳʷ
AMESH,8                    !ͷǵӳʷ
asel,s,,,9                 !ѡֹ֧
aatt,,2                    !
LCCAT,8,11                 !߶
MSHKEY,1                   !趨ӳʷ
AMESH,9                    !ӳʷ

asel,s,,,2,3,1             !ѡֹ֧ڵ
aatt,,2                    !
amesh,2                    !ֹ֧ڵ
amesh,3

Alls                       !ȫѡ
ARSYM,Z,all, , , ,0,0      !Գ
NUMMRG,ALL, , , ,LOW       !ϲ
NUMCMP,ALL                 !ѹ
fini

!*******************************************
/solu
ANTYPE,STATIC              !
outres,,1                  !趨ѡ
NLGEOM,ON                  !η
pred,on                    !ԶԤ
KBC,0                      !Լ
NSUBST,20                  !趨Ӳ
ARCLEN,ON,5                !

lsel,s,loc,y,-hc           !ѡͲ϶
DL,ALL, ,ALL,              !ȫԼ
local,11,1,,,,,90          !ֲϵ
lsel,s,loc,y,60            !ѡṹԳλϵ
lsel,a,loc,y,-60
lsel,u,lcca                !ȥ
dl,all,,symm               !ʩӶԳԼ
alls

lsel,s,radius,,d2/2,(d2+ds1)/4	!ѡ֧
nsll,,1                         !ѡڵ
*get,ntt,node,num,count         !ȡѡڵ
nfyy=fyy/ntt                    !ֱغƽ
f,all,fy,-nfyy                  !ѡڵʩӴֱغ
alls
SOLVE
FINI
!*******************************************
/POST1

SET,list,2                      !бȷٽغ



